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ABSTRACT 



The elastic differential scattering cross sections for 

+ + + + o 

H , + Ar -> H , + Ar were measured for angles of 39 to 

50° with incident lab energy of 100 eV. The apparatus utilized 

the proper ties of axially symmetric, non-uniform magnetic fields 

(large thin solenoid) to focus all scattered ions with vector 

momentum p within Ap onto a detector placed down the magnetic 

axis. Thus, the solid angle was increased by a factor of about 

300 over conventional methods with comparable resolution in AO 

(about 1°). A typical value for (7 (43°) is 8.2 x 10 ^ cm 2 . 

H + 

The H data was compared with data of Magnuson, et.al. , via a 

scaling law and was in excellent agreement. A typical value 

for a (43°) is 2.4 x 10 20 cm 2 . A simple refinement of a 
H 2 

classical model for scattering by Feist relates these cross 
sections by O (0) = P CT (©) ~ Lc 7 (©)] where P is an 

u“*“ _/j^„ TT " r " 



Hr 



ttjTe * 



orientation factor, D the binding energy, and & the internuclear 
separation. This expression relates the cross sections to within 



the accuracy of the experiment thus verifying the classical 
approximations . 



i/ctu y 
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INTRODUCTION 



I_. 

Atomic processes such as the scattering of complex ions 
by heavy atoms are properly described by quantum mechanics. 
Unfortunately, for a many-body problem the mathematics becomes 
untractable. In order to produce useful solutions approximations 
must be introduced. One class of approximations involves the 
use of semi -classical models of the collision process. By 
sacrificing rigor, useful solutions are obtained that still 
provide insight into the physical process involved. 

The classical model investigated here is used to calculate 

+ 

the differential scattering cross section of the ion by a 
spherically symmetric potential. The theory was developed by 
Feist [l] based on the earlier work of Bates, Cook, and 
Smith [2] as follows. First, calculate the differential 
scattering cross section of a single proton from a spherically 
symmetric potential. Then calculate the probability that the 
second proton is oriented with respect to the first such that 
the H* molecule is not disrupted by the collision. The dif- 
ferential scattering cross section for is the product of 
this orientation probability times the cross section for the 
single proton. 

The purpose of the experiment reported here was to test 
the theory as developed by Feist. This was done by measuring 
the differential scattering cross section of H by Argon and 
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H 2 Argon at incident energies of 100 eV and laboratory angles 
from 39 ° to 49 °. Thus, the orientation probability was deter- 
mined experimentally and compared with Feist f s calculation. 
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II. THEORY 



A. THE MODEL 



This classical model is applicable to the scattering of the 

+ 

H^ ion over a range of incident energies (in laboratory 

coordinates) by a spherically symmetric potential such as heavy 

noble gas atoms. The reaction studied was H* + At H* + Ar 

where the underline indicates the fast particle in laboratory 

coordinates. The two protons of the H* ion have an equilibrium 

separation in the ground state of 1.06 angstroms and a disasso- 

ciation energy of 2.648 eV . The lower energy limit where this 

model could no longer be expected to apply would be where the 

incident ion energy becomes comparable with the ion binding 

energy. This requires that E q >> 2.648 eV or E q ^ 3^ eV. The 

upper energy limit where this model could no longer be expected 

to apply is governed by the uncertainty principle. In this 

application it is required that the de Broglie wavelength be 

n h 

much less than the distance of closest approach or r » A. = 

^ o D mv 



lOh 



Now r may be further 
o 



which may be written as r > 

0 ~r^F o 

related to E q by considering a head-on collision with the shielded 
coulomb potential for Argon (this potential is discussed in Section 
II B s and Appendix A ) 0 Considering only the dominant term of the 
potential one can write 



E 

o 




-ar 

o 

Ae 



where A and O' are parameters of the potential. Substituting 



values for the constants one obtains 



E o = 14.6 -§- 
r 

o 

-1.68 r 

and E = 4ll — e ° 

o r 

o 

with E in units of eV and r in units of Bohr radii. The value 
o o 

u 

of E that satisfies these equations is E = 1.13 x 10 eV . This 
o o 

shows that a classical model for the scattering of H* should be 

valid for incident energies up to about 10 KeV. 

. + 

In a scattering event the ion may be disrupted due to the 
following two effects. First, consider the two protons approach- 
ing a spherically symmetric potential with the same impact para- 
meter, b. Using the concept of an impulse approximation, con- 
sider each pr oton- tar get collision separately, yielding a de- 
scription in terms of two two-body collisions. Due to the 
spherical symmetry each proton scattered is confined to a plane 
defined by the scattering axis and a line from the scattering 
center to the proton. The angular difference between the 
scattering planes, A0 will result in a difference in velocities 

of the two protons after scattering such that their relative 

. . + 

energy may be greater than the binding energy and the ion 
will be disrupted. Figure 1 illustrates this process. 

Second, consider the protons oriented such that they will be 
scattered in the same plane but with different impact parameters. 
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Figure I. Protons with same impact parameter but different scattering planes. 



This will result in a difference in scattering angle A 9 and 
therefore again a relative velocity between the two protons. 

The relative energy may again exceed the binding energy so that 
the ion will be disrupted. Figure 2 illustrates this process. 

In the actual case these two processes are not separable 
and disruption may occur due to the combined effect of A0 and 
A9. The necessary condition for the ion remaining intact after 

i 2 

collision is D>%/i(AV) where D is the disassociat ion energy, 

/i = M^/2 the reduced mass of the two protons, and AV is the 

difference in velocities due to both A0 and AO. Whether the 
+ + 

+ Ar — H 2 + Ar reaction will occur is then a function of the 

orientation of the two protons before (or equivalently after) 

the collision. The differential scattering cross section for 

the H* ion, CT (0) may be calculated as follows. For E >> D 
H 2 

treat the collision of the ion as two events consisting of 
the scattering of proton "A" and proton "B" (impulse approxi- 
mation). Calculate the differential cross section for proton 

"A" , CT (0). Calculate the probability, P, that proton "B n is 

H+ 2 
oriented with respect to proton "A" such that D> ^ /!(AV) . It 

then follows that 

a (0) = p a (9). (l) 

h 2 h 

The electron is not treated as a scattered particle in this 
model but may be thought of as simply providing the binding 
ener gy . 



i4 
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Figure 2. Protons with the some scattering plane but different impact parameters. 



B, CALCULATION OF a (Q) 

H + 

The calculation of Cr (0) was done on an IBM 36 O -67 digital 

H + 

computer using a program written by Feist (_ 2 J . The program 

utilizes standard numerical techniques to solve the classical 

differential scattering equations. The potentials used are 

shielded coulomb potentials from Byatt C 3 J of the form 

2 -b r 

V(r) =2L(SCe " ) 

j. n 11 



with SC =1 but C and b otherwise adjustable. Byatt lists 
n n n n 

two potentials for Argon, one with and one without correction 
for polarization. The two potentials yielded cross sections 
for H + that were within 13% of each other for an incident energy 
of 100 eV and scattering angles from 3^° to 50° (laboratory co- 
ordinates). Due to the close agreement of the two potentials, 
only the one corrected for polarization is presented here. 

Byatt indicates the potentials should be accurate to within 5% 

-b r 



for SC e n >.01. For 
n n 

used it was noted that S 

n 

these potentials should 



the region where the potentials were 
-b r 

C e n > .052. Thus it was felt that 
n 

be satisfactory for calculating Cr (9). 

H + 



C. CALCULATION OF THE ORIENTATION PROBABILITY 

+ , ... 

The energy criterion for the ion remaining intact may be 



put as 




( av y 

x max' 




( AV ) 
v max 7 



2 



or 



AV 

max 




( 2 ) 
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Considering the case where the target mass is much greater than 



the mass of the H* ion, the H* ion has energy after collision 
2 

E = M V both in the laboratory and center of mass systems, 
o p o 

Substituting this into equation (2) gives 



AV 



max 



4d 

E 



V 



Now it was shown in Section II A, that AV = f(b, Ab, A0). 
Therefore, for a fixed incident energy and impact parameter, 
there exists a function g such that A0 = g(Ab), This defines 
a closed curve in a plane perpendicular to the scattering axis 
that proton "B ,! must pass through in order for the ion to remain 
intact. The area enclosed by the A0 = g(Ab) curve is designated 
by S . 

The axis of the H* ion is randomly oriented with proton f, B" 
having no preferred position with respect to proton ft A". The 
location of proton "B fl will therefore define a sphere about pro- 
ton "A" of radius £ , the internuclear separation. The orientation 
probability is then given by the ratio of S to one-half the area 
of this sphere 



P = 



2ttj£‘ 



(3) 



The justification for using the area of a hemisphere in the 
denominator of equation (3) is that the H other half" is re- 
dundant. No exact solution has been found for equation (3) and 
the following approximations were made. 



17 



First, it has been assumed that the target is infinitely 

+ 

heavy so that the speed of the ion is unchanged by- the col- 
lision. For the reaction studied experimentally, this is valid 
in that M = 2,032 amu and M = 39*948 amu , For example, a 

100 eV ion scattered at 45° has an initial speed of 9*80 x 10 

m/sec and a final speed of 9-75 x 10^ m/sec. It is further 
assumed that the two protons will both have this same speed after 
collision and that disruption will occur due to difference in 
direction only. Feist L 1 J showed that the disruptive effect 
of the difference in the magnitudes of velocities of the two 
protons is much less than the disruptive effect of their dif- 
ference in direction. For the case where A0 = 0, and assuming 
A0 is small, AV = V [ A0 J . w ~ . Using equation (2) this 
becomes _ . 



^^max^A0 = 0 



4d 

E 

o J 



(4) 



Figure 3 a illustrates this case. For the case where A0 = 0 
(or equivalently Ab = 0) and assuming A0 is small, 

AV = sin 0 [A0_^_] ^ _ Q . Using equation (2) this 



max o 
becomes 



max' 



U0 maAb = o = [r f ii^r 



(5) 



Figure 3b illustrates this case. Using the classical definition 



for differential scattering |2rrbdb! - | CT (0) 2 tt sin 0 d0| , it 

H 



is seen that 



Ab 



a ( 0 ) s in 0 
H + 



max 1 A0 = 0 



[ao ] 

I max J 



max J A0 = 0 5 



18 




Figure 3. a. Proton velocity difference due to difference in impact 
parameter. 






* V Q sin 8 




[ A ^ max ] 



Ab s O 

Scattering axis into paper 



Figure 3. b. Proton velocity difference due to difference in 
scattering plane. 
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which, upon combining with equation (4) becomes 



Ab 



* +( Q ) 

H 



sin © 



max A0 = 0 



4d 

E 



( 6 ) 



Equations (5) and (6) give maximum values for either Ab or A0 
+ 

such that the H 0 ion is not disrupted. The A0 = g(Ab) curve 
must therefore be enclosed by these maximum values. An upper 
bound on the area, S," enclosed by the A0 = g(Ab) curve is there- 
fore 



Subs t itut ing 
equation (3) 
found to be 



and 



S = [2(Ab ). rt _ J[2b(A0 ).. J. 

max max'A0 - 0 v max'Ab = 0 

equations (5) and (6) into the above and using 
an upper bound on the orientation probability is 



P 

max 





(7) 




max 



8d 

rrJl 2 E 



[° + ( 6 ) 3 2 - 

H 



( 8 ) 



Due to the simple form of equation (8), it is calculated numeri- 
cally in the same computer program used to calculate 0“ (©). It 

H + 

should be further noted that D is the disas socia tion energy for 

ground state and JL the corresponding equilibrium distance. Ex- 

. + 
cited electronic states for the molecule are either unbound 

or have very shallow minima in their potential curves. Therefore, 

if the ion is excited to higher electronic states during the 

collision process it is almost certain to be disrupted due to 
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velocity differences after the collision. If the ion is 

excited to higher vibrational levels, then D would be decreased 

and l increased. The molecule could still remain intact but 

would have a smaller value of P and hence O’ (9). 

H 2 

D. COMPUTER CALCULATION OF O (9) AND a (0) 

H + 

Figure 4 displays the results of the computer calculations 
for the differential scattering cross sections. One set of 
curves is for the + Ar + Ar reaction at incident energies 

of 100,300, and 500 electron-volts. The second set of curves 

+ + . , 
is for the + Ar H Q + Ar reaction also at incident energies 

of 100, 300, and 500 electron-volts. The curves were com- . 
puted using equation (8) and represent a predicted upper limit 
on the scattering cross section. 

It is noted that all of the curves are monotonical ly de- 
creasing functions of scattering angle and reflect no atomic 
structur e . 
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Figure 4. Log of scattering cross section vs. scattering angle for 
H* + Ar-» H * + Ar and H + + Ar-*H + + Ar. 
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III. EXPERIMENTAL APPARATUS 



A. GENERAL DESCRIPTION. 

The experimental apparatus used to measure the cross sections 

+ + + + ... 
for the H , H 2 + Ar H , + Ar reactions is discussed in detail 

by Bush L 4] - An overall schematic of the apparatus is shown in 

Fig. 5 • The operation of the apparatus was as follows. Three 

*4* -f* *4* 

species of hydrogen ions (H , H 9 , H ) were produced in the duo- 

2 3 

plasmatr on ion source and were accelerated electrostatically 
to the desired energy. The ions entered a mass analyzer where 
the desired ion species was bent through a 45 angle to travel 
down the geometric axis of the apparatus. The selected beam was 
focused and entered the scattering cell where the target gas, 
in this case Argon, was contained. The scattering cell was 
centered in a large focusing magnet that produced a known non- 
uniform but axially symmetric magnetic field. The trajectory of 
an ion with a given vector momentum could then be calculated and 
the detector was appropriately placed to collect this ion. Due 
to the finite detector size, ions scattered within a finite AO 
(on the order of 1°) were also detected. It is further noted 
that essentially all the ions scattered at an angle 0 were col- 
lected regardless of their azimuthal angle of departure. Cross 
section measurements at different scattering angles could be made 
in two ways. Either the detector could be fixed on the magnetic/ 
geometric axis and the focusing magnet current varied to select 
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FIGURE 5. SCHEMATIC OF THE EXPERIMENTAL APPARATUS. 



the desired scattering angles, or the focusing magnet current 
could be fixed and the detector moved along the magnetic/geometric 
axis to select the desired scattering angles. In this experiment, 
the latter method was used. The detector system consisted of a 
magnetic electron multiplier whose output was fed into an elec- 
trometer. The entire system was maintained at a vacuum on the 
order of 10 ^ torr during operation by two diffusion pumps. 

B. THE DUOPLASMATRON ION SOURCE. 

The duoplasma tr on ion source was constructed and described 
in detail by Carter L 5 J - Figures 6 and 7 display the construction 
of the source and the circuit diagram for its operation. The 
operation of the source basically consisted of maintaining the 
plasma discharge in the Z-electrode cavity and extracting ions 
at the desired energy by placing this potential on the anode. 
Focusing the extracted beam was accomplished "by an electro-magnet 
surrounding the Z-electrode cavity and a conventional Einzel lens 
and deflector stack. 

To increase the beam intensity an acceleration-deceleration 
system was utilized. This consisted of extending the extraction 
lens, , to enclose the entire focusing stack. This enclosure 
was called the "sock". (See Fig. 7 «) By floating the sock at a 
negative potential V ( typ ical ly -500 volts) and by increasing the 
potentials of the interior lenses and deflectors by the same 
amount, the extracted ions transversed this region with an energy 
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FIGURE 7. 

CIRCUIT DIAGRAM FOR DUOPL ASMATRON 
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1 

4 

(E + V)e, yielding a higher space charge limit for the beam 

density , Upon exiting from the sock, the hydrogen ions were 

decelerated to energy Ee. This arrangement increased the beam 

intensity at the scattering cell by a factor of about 40 . 

A typical maximum beam intensity for mass analyzed protons 

of 100 eV delivered to the scattering cell with 0 o 5 cm entrance 

“7 

aperture at 80 cm from the source was 1 x 10 amperes . For 

+ ... 

ions under the same conditions the maximum beam intensity was 

-7 

on the order of 2 x 10 amperes. The energy spread of the beam 
was estimated to be about 2 eV for a 100 eV beam. 

C, THE MASS SPECTROMETER. 

The mass spectrometer used to select the desired ion was 
constructed, calibrated and described in detail by Strohsahl L6] . 

Figure 8 is a schematic of its arrangement in the overall appa- 
ratus. The mass spectrometer utilized an electromagnet to pro- 
duce a magnetic field of up to 4 kilogauss confined between two 
40° wedge-shaped pole pieces. This was sufficient to mass 
analyze atomic masses up to 40 with energies up to 1 KeV. The 

resolution was about 40 . A bending angle of 45 ° was 

Am 

utilized and this altered the selected ion beam trajectory from 
the axis of the duoplasmatr on to the geometric axis of the main 
scattering apparatus. Operation of the mass analyzer consisted 
of merely of setting the correct magnet current to focus the 
desired ion into the scattering cell. The pole faces and region 
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Figure 8. Schematic of Mass Spectrometer. 
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where undesired species were "dumped™ were coated with aquadag 
to prevent static charges from building up. Even so, this region 
required thorough cleaning after about 20 hours of duoplasmatr on - 
analyzer operations. 

D. THE SCATTERING CELL. 

The scattering cell was designed by Bush L 4] . Figure 9 shows 
a schematic of its construction. The functions of the scattering 
cell are to (1) allow the collimated incident ion beam to pass 
through the scattering region and be measured by the beam collector, 
(2) contain the target gas in the scattering region and measure 
its pressure* and (3) allow scattered ions to escape over as large 
a range of scattering angle and azimuthal angle as possible. The 
design of the scattering cell was a compromise between these re- 
quir emen ts . 

The incident beam enters the scattering cell through an 
adjustable aperture (diameter from 0.1 to 1.0 cm) which was 
usually set at 0.5 cm. The beam passed through the front section 
which contains the scattering region. Now the great majority of 
the incident ions will not be scattered at all and will continue 
on into the rear section where they are collected either by the 
rear section wall or the beam collector. The incident beam was 
measured by two Keithly Model 4l0 Micro-Microammeters. The beam 
collector, which collects about 90 % of the beam, had three grids 
that prevented slow ions from being collected and suppressed 
secondary electron emission. The beam collector assembly could 

be swung off axis allowing a clear path through the scattering cell 

for system alignment and detector calibration. 
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Figure 9. The Scattering Cell and Beam Collector 



The target gas was bled into the scattering region through 

a Granville-Phillips variable leak valve. The pressure was 

measured by a VG-lA ion gauge manufactured by Vacuum Industries. 

Since operation of the ion gauge interfered with the incident 

beam measurement , the pressure was measured before and after 

each run and averaged. Vacuum Industries lists the correction 

factor for Argon (relative to dry air) as 0.617. That is, 

P . = o 0617 X P. . 

actual indicated 

The rear section of the scattering cell was adjustable along 
the axis of the apparatus and thus regulated the target thickness 
between 0.1 and 0,8 cm. This was set at 0.3 cm for the present 
experiment. 

The scattering angle limits were defined by the slope of the 
front and rear section shoulders of 49° and 36 ° respectively. 

This range of scattering angles could be extended somewhat by 
considering an effective target thickness. Bush L4*J describes 
this process. Ideally, ions scattered into all 360 ° of azimuthal 
angle could be collected at the detector. However, interference 
by mechanical support required for the scattering cell reduced 
the available azimuthal angle to 290 °, 

E. THE FOCUSING MAGNET. 

The focusing magnet consisted of five concentric spools of 
copper wire in a par allel-ser ies arrangement. Each spool was 
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water-cooled by a ^ inch copper tubing jacket. A Sorenson 
Nobatron DCR 150-15 power supply maintained the desired current 
setting for periods of several hours with essentially no drift. 

The magnetic field on the center axis was carefully surveyed 
by Kelly L 7 J using a Hall probe. Figure 10 shows the result of 
these measurements. Kelly also showed that the magnetic field 
was symmetric about the coil axis to within 0.1%„ This allowed 
the axis of the apparatus to be aligned with the coil axis so 
that the magnetic axis coincided with the geometric apparatus 
axis . 

In order to better understand the function of the focusing 
magnet, it is instructive to consider the simpler case of a 
uniform axially symmetric magnetic field. This would correspond 
to the entire apparatus being enclosed in a long solenoid. For 
an ion of charge q with momentum p and scattered at an angle 0 
the trajectory is a helix such that after one orbital revolution 
the ion will again intercept the axis at some distance Z q . 

Using these boundary conditions the Lorenz force equation 

» 2 tt 

F = q v x B can be integrated to give Z = — p cos Q. A de- 

o qB ^ 

tector placed at Z q will then detect the scattered ion. (See 
Fig. 11.) This equation also shows the two ways of detecting 
scattered ions over a range of 0 and p. Either maintain B con- 
stant and vary the detector position or keep the detector fixed 
and vary B. 
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Figure 10, Magnetic Field Strength at Various Distances 
Along the Magnetic Axis for Various Currents 
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Figure 11. Orbit of Charged Particle in Uniform Magnetic Field 



Now in order to measure a differential scattering cross 
section, the solid angle that the finite size detector "sees" 
must be known. For the axially symmetric case this is 
d Cl = 2 tt sin 0 d 0 . To evaluate d© consider the various para- 
meters affecting Z q , the axis intercept distance, A functional 
relationship exists such that; 



wher e 



Z q = f(E,0,B,/3,7?,C) 

E = energy of incident ion 



0 = angle of scatter 

B = magnetic field 

p = radial target displacement from the system 
axis 

Tj = displacement of target atom from the 

scattering center along the system axis 

C = constant for each reaction including mass 
ratios, inelastic energy losses, etc. 

Differentiating this expression gives 




For a specific measurement AB and AC are zero. Since analytic 
expressions for the various derivatives are not known, they must 
be evaluated numerically by varying E, 0, p, and 7} . Bush, Cook, 
Heinz, and Rodeback [ 8 ] estimate that Ip and tor) which are due 
to the finite volume of the scattering region affect A0 by no 
more than 10%. Approximating A p and A 77 as zero the above ex- 



/d© ' 

pression can be rearranged to A© = / - 



- (If) 



AE 



36 



The effect of a small ( 2 %) energy spread was studied and it was 
concluded that the "effective" A© was not changed. The angular 
resolution of the apparatus was degraded somewhat in that not 
all particles scattered within A© of © were collected but this 
was compensated by the fact that some particles scattered outside 
of A© were collected. The solid angle can therefore be approxi- 
mated as ✓ \ 

dfi - 2n sin 0 j ( ) 1 AZ o 

AZ o 

where _+ — is the maximum deviation of the trajectory inter- 

cept from Z q that will still allow the particle to be collected 
by the detector. AZ q can be calculated by considering the de- 
tector geometry and the angle © o with which the particle inter- 
cepts the Z-axis. As an example of how the approximate expression 
for the solid angle is applied if an analytic expression for 

Z q = f(©) is known, consider again the uniform axially symmetric 

2tt 

magnetic freld. Recalling that Z q = cos ©, it follows that 

dCl ~ - 2*5 az . 

- p o 

For the non-uniform axially symmetric field, the procedure 
for computing trajectories and solid angles was essentially the 
same as for the uniform case except that the calculations were 
done numerically on an IBM 36O-67 digital computer. The program 
that accomplished this is explained in Appendix A. A very brief 
outline of the program is as follows. 
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(1) Express the magnetic field along the axis as a 12th 



I is the magnet current. 

(2) Using this polynomial, a vector potential is found. 

(3) Combining the vector potential with the Lorentz force 

equation gives a 2nd order non-linear differential 

equation. This is expressed as two simultaneous, first- 

order differential equations which are integrated by 

the DHPGC Subroutine of the Scientific Subroutine 

Package for the IBM 36 O -67 computer. The trajectory, 

Z , and Q are now known, 
o o 

(4) The above procedure is done over a selected range of 
scattering angles and a polynomial expressing 0 = f(Z Q ) 

is generated. 

(5) From 0^ and the detector geometry, AZ q is calculated. 

(6) The solid angle can now be calculated as 



is present due to the construction of the scattering 
cell o 

Figure 12 shows a series of trajectories computed for H 
ions with incident energy of 100 eV elastically scattered by 
Ar gon . 



^ n 1 

order polinomial. This gives B = (I) £ C Z where 

z v 1 1 n 
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focusing magnetic field with I = 4.8 amps. 



F . THE DETECTOR . 



The detector system employed in this experiment consisted of 
a Bendix Model 306 Magnetic Electron Multiplier (MEM) which was 
fed into a Keithly Model 640 Electrometer t Figures 13, l4 and 
15 show the detector geometry, a schematic of operation, and 
a circuit diagram. The principle of operation of the multiplier 
was as follows. An incident ion was accelerated to about 1500 eV 
and ejected an electron from the tungsten cathode upon impact. 

A crossed electric and magnetic field (alnico permanent magnets) 
imparted a cyclical trajectory to the electron(s) down the dynode 
strip. The glass dynode and field strips had a high resistance 
coating of tin oxide and antimony that sustained the electric 
field. Each time the electrons struck the dynode strip multi- 
plications occurred so that when the electrons were collected by 

4 

the anode, a gain of about 10 had been realized. 

Unfortunately the MEM had to be operated in the focusing 
magnetic field. This resulted in a sideways drift of the elec- 
trons moving down the dynode strip so that all of the electrons 
generated were not collected at the anode. This in fact made 
the MEM gain a very sensitive function of its position in the 
focusing magnetic field. The procedure used was to calibrate 
the MEM under the same environment that existed during the actual 
cross section measurements. Specifically, this consisted of: 

( 1 ) measuring the gain using ions with a focusing magnet cur - 
rent of 6.6 amperes with the detector position, Z q , varied from 
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Figure 13. Dector Aperature Geometry 
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FSI 



Field 

Strip 



FSO 



Input differential 
(250 Volts) 




Differential 
(250 Volts) 



Dynode Voltage 
(1500 Volts) 



Figure 15. 



Circuit Diagram for Electron Multiplier. 



43 



30.5 to 47«5 centimeters from the scattering cell, (2) measuring 
the gain using H ions with a focusing magnet current of 4-8 
amperes with the detector position, Z q , varied from 28.5 to 42.5 
centimeters from the scattering cell. Since the focusing magnetic 
field was a known function of current and position, the two gains 
are plotted against magnetic field in Fig. 16 . The gain for H+ 
ions was measured before and after the cross section measurements 
and was found to remain constant. The gain for H + ions was 
measured only once. 

Due to the steep slopes of the gain curves, a small error 
in setting detector position or focusing magnet current will 
result in a large error in the detector current reading. This, 
coupled with the tedious task of calibrating the detector for 
every species, energy, and magnetic current that a cross section 
measurement is to be made, made this detector system less than 
ideal. For further experimental work with this apparatus, a 
different detector system is being considered. This may consist 
of counting the output pulses from the MEM rather than integrating 
them with an electrometer. 
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Figure 16. Multiplier gain vs. magnetic field. 
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IV. EXPERIMENTAL RESULTS 



A. LABORATORY CROSS SECTION MEASUREMENT. 

The equation used for calculating a differential cross section 
in the laboratory was 



V G 



^ ^ ^ I x p sc x Corr .Factor )x t x d £2 ( 9 ) 



2 

cm . 



wher e 



Ip = detector current measured by an electrometer 



G = MEM gain 

I. = Incident ion beam. Sum of the currents to 
1 

Beam Collector and Rear of Scattering cell. 

1 6 3 

p = 3 .54 x 10 par ticles/cm -tor r based on ideal 
gas law. 

P = Pressure in scattering cell as measured by 
sc 

VG-1A ion gauge 

Corr . Factor = .617; corrects ion gauge reading for Argon 
t = target thickness. Normally equal to . 3 cm 
in this experiment 

d Cl ( 0 ) = Solid angle computed by digital computer. 

A sample calculation for the H + Ar reaction showing typical 
experimental values is 



a (Ll 9°) - ( 1 - 8 5 x l(T 9 amp.)/1.27 x 10 ^ 
ct ) - _» 16 -3 



H 



Q 

(1.4l * 10 amp 0 )(3.54 x 10 x ^cm J -torr x ) 



X 



( 1.38 X 10" 3 torr X .6l7)(0.3 cm)(0.109 sterad) 
-17 2 

= 1.07 x 10 ' cm 
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B. MEASUREMENT OF THE H + + Ar - H + + Ar CROSS SECTION. 

The differential scattering cross section for this reaction 
was measured from 39-2° to 50°. Five runs were made and the 
data points are displayed in Fig. 17 along with the computed 
cross section using the Byatt potential. Figure l8 shows the 
average of the 5 runs along with the computed cross section. 
Important parameters for the measurements included: 

(1) Incident energy of H + ions equal to 100 eV. 

(2) Focusing magnet current equal to 4.8 amps. 

(3) Solid angle typically equal to 0.1 steradians. 

(4) Argon pressure in scattering cell typically equal to 

-3 

1 x 10 torr . 

Figure l8 shows that the experimental curve, while having 
about the same overall slope as the computed curves , differs in 
magnitude by a factor of about 6. This difference, while not 
catastrophic, was greater than expected. First, there may have 
been a systematic error in the experimental apparatus. A pos- 
sible systematic error could have been an error in measuring the 
MEM gain. The gain was measured with an incident beam several 
orders of magnitude greater than the beam of scattered particles. 
Thus, if the MEM gain were a (slowly) increasing function of input 
current, the measured gain would be greater than the gain for the 
smaller beam of scattered particles. This would result in the 
measured cross sections being too small. 
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Figure 17. Experimental measurements of _H* 4 Ar -► _H_ + + Ar at 
100 eV. 
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H*+ Ar — ► H + + Ar at lOOeV. 
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Second, the potential may have been inadequate for this 
application . Varying the potential parameters revealed that the 
computed cross sections were relatively insensitive to small 
changes in the potential parameters. For example, varying the 
dominate terms in the potential by 30 % such that the collision 
was made softer decreased the computed cross sections by a factor 
of only about 2. It was felt that the potential parameters were 
at least this accurate. It is noted, however, that the potential 
is Hartree-like and therefore will reflect no atomic structure 
or inelastic process. At the large angles where this experiment 
was conducted, the collisions are relatively hard so that the 
Argon atom may well be excited. The H + ion would then have 
significantly less energy and momentum after the collision com- 
pared to the elastic case. Since the focusing magnetic field 
acts as a- momentum analyzer, only elastically scattered H + ions 
were collected by the detector. Thus, the presence of an 
inelastic process could decrease the elastic cross section by a 
factor of 6 and not be accounted for in the Hartree-like potential. 

The work of Aberth and Lorents [ll] tends to support the 

second analysis. Aberth and Lorents measured the scattered signal 

profile of He on Ar at 10 and 400 eV with an energy analyzer. 

+ 

The result shows that for He on Ar , an inelastic event is more 
likely than elastic scattering by a factor of about 2%. The energy 
lost by He + in inelastic events was on the order of 10 to 15 eV. 
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Applying this to the case of H on Ar , it seems reasonable to 

expect protons to be a more efficient projectile for exciting 

+ - 
Ar than He and therefore could account for the factor of 6. 

Furthermore, inelas tically scatter protons with energy after 

collision of 85 - 90 eV would not be collected by the detector. 

The choice between the computed and experimental cross 
sections was accomplished by comparing the data in Fig. l8 with 
data supplied by Magnuson, Carlston, and Snyder [ 9 ]* Magnuson 
et al also measured the cross section of H + on Ar although their 
measurements were made for scattering angles less than 20° and 
energies ranging from 150 eV to 800 eV . Meaningful comparison 
of the data was possible by using a reduced plot outlined by 
Smith, Marchi, Aberth, Lorents, and Heinz [l0]. This reduced 
plot consists of plotting T = EO versus Q = Osin 0 0(9). By 
the use of this scaling law, cross section curves obtained at 
various energies are reduced to a single line. The reason is 
that the product EO represents essentially the same encounter 
with the atomic potential. By comparing the data in this manner 
it was found that the two experimental curves agreed to within 
a few percent and the unaveraged data infact overlapped. Thus, 
it was concluded that the experimentally determined cross sections 
were probably correct. It is further noted that the data of 
Magnuson, et_.al_. covered a much larger range of T = EO than this 
experiment and does reflect atomic structure such as discussed 
in the preceeding two paragraphs. 
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Co MEASUREMENT OF THE H* + Ar + Ar CROSS SECTION. 

The differential scattering cross section for this reaction 
was measured from 38.8° to 49«7°. Five runs were made and the 
data points are displayed in Fig. 19. The average of the 5 runs 
is displayed in Fig*, 20. Important parameters for this measure- 
ment included; 

(1) Incident energy of ion beam equal to 100 eV. 

(2) Focusing magnet current equal to 6.6 amperes. 

(3) Solid Angle typically equal to 0 o l steradians. 

(4) Argon pressure in the scattering cell typically equal to 

-3 

1 x 10 torr . 

From Fig. 20 it is seen that the overall slope of the experi- 
mental curve is essentially flat. That the cross section does 
not decrease significantly with increasing angle is possibly due 
to the following instrumental effect. Some particles scattered 
outside of the computed solid angle may reach the detector through 
multiple collisions with the walls of the apparatus and still have 
sufficient energy to overcome the detector grid bias. The sig^ 
nal due to these extraneous particles could not be subtracted 
out in that they essentially behaved like the "real" particles 

scattered into d £7. That is, the signal due to extraneous 

. . . . + . 
particles was proportional to the incident beam intensity and 

to the pressure in the scattering cell. Normally, this extraneous 
signal is negligable compared to the signal from particles scat- 
tered within d£lo However, in this case at the larger angles the 
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4* CM 



detector signal was quite small, typically on the order of 
-12 

0.3 x 10 amperes. Therefore, perhaps up to 50 percent of 
the measured signal could have been due to these extraneous 
particles. This would also account for the wider dispersion of 
data at the larger angles. 

D. EXPERIMENTAL RESULTS AND THE CLASSICAL MODEL. 

Since it could not be assertained with certaintly as to 
whether the computed or experimental cross section for H + on 
Argon was more accurate, the analysis of the classical model 
was done from both points of view. 

First, it was assumed that the computed cross section for H + 
was correct. The difference between the experimental and com- 
puted cross section was then attributed to a systematic error 

+ 

and the experimental H data was normalized to fit the computed 
curves. This was done by moving the H + data up by a factor of 6. 
Now the on Argon cross section measurement had the following 
similarities with the H measurement: 

(1) Approximately the same focusing magnetic field (see Fig. l6 ) . 

(2) Approximately the same excursion in detector position. 

( 3 ) Approximately the same calculated solid angle. 

( 4 ) The gain curves for both measurements have about the same 
shape (see Fig. 16 ) . 

It was therefore concluded that if a systematic error were pre- 
sent, it should have affected the measurement by about the 
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same magnitude. Therefore, the data was also increased by a 

factor of 6 . The results of this normalization are shown in 

+ 

Fig. 21. It is seen that the on Argon experimental data pre- 
sented in this manner lies below the predicted curve by a factor 
of about l4« 

Second, it was assumed that the experimentally measured H + 

on Argon cross section was more accurate than the computed. 

Since the classical model predicts that the cross section for 

H* is proportional to the square of the H + cross section (Eq. 8), 

+ 

the predicted curve for on Argon m Fig. 21 would be moved 
down by a factor of (6) or 36. Figure 22 illustrates this. In 
this case, it is seen that the experimental measurement for H+ 
on Argon lies below the cross section predicted by the classical 
model by a factor of about 2 to 3. 
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Figure 21. Analysis of experimental measurement assumming a 
systematic error. 



57 




Figure 22. Analysis of experimental measurement assumming no 
systematic error. 
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V. CONCLUSIONS 



The experimental results given above can be analyzed in two 
distinct ways: (1) Assume that the theoretical H + + Ar curve is 
correct and that a systematic error was present in the experi- 
mental measurements. ( 2 ) Assume that the experimental H + + Ar 
cross section is correct. We present below the conclusions based 
on both assumptions and our reason for believing that method (2) 
is more likely to be correct. 

If a systematic error were present and affected the measure- 
ments of both the H and cross sections by the same magnitude, 

then it was seen that the classical model predicts a cross section 

+ . 
for too large by a factor of about lq . That the actual cross 

section for is so much smaller than predicted is probably due 

to the two effects previously discussed. That is: (1) the 

orientation probability wa-s computed by considering the maximum 

area enclosed by the A0 = g(Ab) curve and ( 2 ) the model neglected 

the possibility that the molecule might be excited during the 

collision process. That the predicted cross section was so much 

greater than observed would perhaps make the classical model 

useful only for predicting upper limits of cross sections . 

However in view of the agreement of our H on Argon data with 

Magnuson et.al. [ 9*1 , this is considered the less likely of the 

two conclusions. 

If both experimental measurements are correct, then it is seen 
that the classical model predicts an H+ cross section too large 
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by a factor of about 2 to 3* This would make the classical model 



useful for predicting H 0 cross sections on an absolute basis . 
This is considered to be the more likely of the two conclusions. 
Assuming that the classical model is essentially correct the 
following simple refinement is proposed: Recalling equation (3) 

it was seen that the orientation probability could be written as 

where S was the area enclosed by the A0 = g(Ab) curve. Since S 
could not be determined exactly, an upper limit was obtained by 
considering the maximum values for A0 and Ab respectively when 
the other variable was zero. This in fact gave four points on 
the A0 = g(Ab) curve, namely: 

A0 = 0, Ab = + Ab 

max 

A0 = 0, Ab = - Ab 

max 

Ab = 0, A0 = + A0 

? max 

Ab = 0, A0 = - A0 

max 

A first appr oximat ioh “ to the A0 = g(Ab) curve can therefore be 
obtained by simply connecting these points with straight lines. 
The following sketch illustrates this. 
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The area enclosed by this first approximation is 

S = %[2(Ab )A0 = o] [ 2b( A0 )Ab = o] = S 

v max 7 v max 7 max 

which leads to P = ^ P 

max 

and a + (9) = [a + (0)] 2 . (9) 

H 0 TT^ E H 

^ O 

Equation (9) predicts the experimentally observed cross section 
+ 

for to within the accuracy of the experiment. 
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APPENDIX A 



1. Computer Program for (T + (0) and a (9) 0 

H H2 



This computer program solves the following set of equations 
for classical differential scattering in the center of mass 
coordinate system. 



* + ( 0 ) 



b | db i 
Sin9 d9 



wher e 



Q = rr - 2bU 



U = 




r 

o 



largest positive root of 1 



V(r ) b_ . 

i 2 2 

%UV r 



0 



with jl = reduced mass of the system 

V q = velocity of incident particle in the laboratory 
system 

V(r) = scattering potential 

In this case the scattering potential as given by Byatt (7 3^1 

1 C^R 8R VR 

for Argon converted to electron volts is V(R) = “(Ae + Be H +C ) 

K 

with R in units of Bohr radii (0 o 529 2%) where 



a = 411.4 


a - 


- 1.676 


b =-117.5 


P = 


- 3-127 


c = 195.9 


V = 


- 9.625 
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An outline of this program is as follows: 

1. Read in the desired incident energies, mass of proton, mass 
of target atom, disassocia t ion energy, and potential parameters. 

2. Select an impact parameter, b. From an initial value, b is 

increased in units of .05 a where a is the first Bohr radius. 

o o 

3. Solve for the largest positive root of 



v ( r D ) b 2 

1 - - ^ = 0 using the Newton-Raphson method. 

/IV / 2 r 

o o 



4. Evaluate the integral 



- u - j 



vtil 



^ V o 



using the Gauss-Mahler quadrature formula 




dr 



U 




H> 9<*1 



) + ( COS 



§)g(x 2 



) +( COS 




) 



where g(x^) 



1 



V ( r /x.) 
x o i' 




b 



2 







-h 



r r 3 tt 5^ 

with x = cos — , x = cos — , x = — . This procedure was taken 

_L _L nS ^ JL j _L 

from Bates, Cook, and Smith [ 2 ]. 



5. From 0 = tt - 2bU, evaluate | TT I = |2b + 2b is 



dU 



dU . 



db 



db 



db 

obtained from the analytic expressions 

V( r o ) ,2 

U = f(r ,b) and 1 - — - — ~ = 0 in a straight-forward 

/i/2 V 2 r 
“ o o 

dr 

manner using the chain-rule of differentiating — = — + 



6. Solve for a (0) 
H + 



_b j db 

s inQ ' d0 
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8d 



7« Solve for the orientation probability, P = — 7 ? — CT (9) 

n£ 2 E H + 
o 

8. Solve for or (9) = P a (9) . 

Hi H + 

2 



9. Write out incident energy, b, r , 9, (7 (0), P, or (9). 

H + HJ 

Computations are done in units of Bohr radii, but CT (9) and 

H 2 

2 

CT (9) are converted to cm at the end of the program. The 
H 2 

following is a list of key symbols used in the computer program 

DENER = Dfsassocia ticn Energy , D, in units of eV . 
RHO(J) = Impact parameter, b, in units of Bohr radii. 
EV(K) = Incident ion energy in units of eV . 

definition 



ET = 



M 



H 



V ( R ) R 

V = — convenient form of the potential in 

this program 
dV 

DV = 

DEG(J) = scattering angle in degrees converted to 

, laboratory system 

SIGMA = Or (0) 

H 



cs = a (0) 

H 2 
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2 . Calculation of Ion Trajectory and Solid Angle » 

As discussed in Section III E, the scattered ion trajectories 
and solid angles must be computed numerically,. Since much of the 
computer program that accomplished this is peculiar to the Naval 
Postgraduate School IBM 36 O -67 computer installation, only an 
outline of the program is presented here. 

Recall that the axial component of the magnetic field on the 
magnetic field axis was measured and expressed as a polynomial 

of the form 13 

B ('z , r - 0) = I E C z n 
z v . ' 1 n 

Bush [4J shows that a vector potential for this axially symmetric, 
non-uniform magnetic field can be written using this expression 
for B^(z,r = 0) as s 



A (r z) = e 

V ’ ' ■ n=0 n!(n+l)j 



(|) 



2n + l & 



2n 



B (z 5 r = 0) 



dz 



2n 



A r (r ,z) = 0 
A z (r ,z) = 0 

Using this vector potential in the Lorentz force equation 

—4 — » — ♦ — » — > — f 

F = qvx B=qv x(V x A) , 



it follows that: 




fd0 a(t 

dt 3r 


> 




<i) 




/ dz 


\/ aA «l 


| 1 dr 3 < r V 


= - q 


(dt 


A dz / 


r dz dr 


= q r 


d0 

dt dz 


(3) 



( 2 ) 
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. Us ing this , 



Equation (2) can be reduced to — - = - q — 

1 v ' dt ^ mr 

equations (1) and (3) can be combined to eliminate the time 
dependence to 




This is a second-order, non-linear differential equation whose 
solution describes the trajectory of the scattered ion in the 
magnetic field. In order to integrate it numerically, it is 
expressed as two first-order, simultaneous differential equations 
as follows. 

Let y x = 



r 




dz dz 



dA 



0 



f. = 



0 dz 



— - A 2 

2 A 0 

q 



SA 



0 



f = 



_ 0 dr 



V 



Then equation 



(4) may be equivalently written as 
dy i 





f 2 V 2 + f i y 2 



( 5 ) 

( 6 ) 
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Equations ( 5 ) and ( 6 ) were integrated by the DHPGC subroutine 
which employs a Hamming v s modified predictor -corrector method. 
The solution to ( 5 ) and ( 6 ) yields the trajectory and hence Z 
and 0 . Now in a single computer run, the trajectory equations 
were solved a number of times for various scattering angles. 
Typically this would be for 0 from 35 ° to 55 ° in 1° increments. 
A second subroutine 'was then called on to generate a seventh 
order least squares polynomial fit of 0 = f(Z Q )o From 0 q and 
the detector geometry, AZ q was computed. The solid angle could 
then be computed by the expression 



aa ~- si " 0 Iff- I 4Z o • 



The computer output then tabulated 0, Z q? AZ q , 
well as a plot of the trajectory (see Fig 0 12) 



0 , and d Cl as 
o 

if desired. 
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